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ABSTRACT 



Galaxy cluster gas mass fraction vers us redshift data (jAllen et al.ll2004l ) and 
Hubble parameter versus redshift data (jSimon et al.l 120051 ) are used to jointly 
constrain dark energy models. These constraints favor the Einstein cosmological 
constant limit of dark energy but do not strongly rule out slowly-evolving dark 
energy. 



Subject headings: cosmology: cosmological parameters — cosmology: observa- 
tions — X-rays: galaxies 



1. Introduction 



Type la supernova apparent magnitude versu s redshift data now favor nonzero dark en 



Clocchiatti et al. 


2006; 


Astier et al. 


2006 



Riess et al.l 120071 : iDavis et al.l 120071 ). Consistent with this, cosmic microwa ye background 



anisotropy measurement s indicate that the Universe is spatial 



200ll : iDurrer et all 1200.1 IPage et al.l 120031: ISoerge 



assumed to be a constant, see, e.g., I Wright 



2006 



y flat (see, e.g.. iPodariu et al. 



et al.ll2007l, if the dark en ergy density is 



Wang fc Mukherjed 120070, and, in con- 



junction with the low observed nonrelativistic matter density (see, e.g.. lChen fc Ratrall2003bl : 
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Spergel et al.ll2007l ). imply that dark energy (DE) accounts for 
ergy budget. 



70 % of the Universe's en- 



A number of explan ations have b een proposed for the DE phenomena. DE might be 



a cosmol ogical constant ([Peebles 1 19841 ) or it could be a dynamic scalar fi eld with negative 
pressu re ( Peebles fc Ratraf 1988 )_ I For recent dark en ergy reviews see, e . g^ Peebles fc Ratra 
j2003UPadmanabhanl (boost ). ICopeland et al.l tood ). and iNobbenbfil (hood ). 



Since different DE models make different predictions for the expansion history of the Uni- 
verse and for the growth of perturbations, DE model parameters can be constrained by using 
available cosmological observations. O bservations such as Type la supernova (SNIa) appar- 



ent luminosity vers us redshift (see, e.g. JNesseris fc Perivolaropouloal200a ; Ijassal et al 



2003 



Barger et al.l 2007 ); cosm ic microwave background anisotropy (see, e.g., Mukherjee et al. 



2006 



Spergel et aD l2Q07h: the angular size versus redshi f t relation for quasars an d radio 



Chen & Ratrall2003al : IPodariu et al.ll2003l : iDalv fc Diorgovskill2006r): strong 



sources (see, e.g.. 

gravitational lensing by a foreground galaxy or clu ster of galaxies (see, e.g., IChae et al. 



2004; Alcaniz et al. 



surements (see, e.g.. ISeljak et al. 



2005 



2005; Fe deli fc Bartelmannl|2007f); and various large-scale stru cture mea- 



Tegmark et al 



200 



6; 



Eisenstein et al 



'ercival et a 



20051 : IWanefeood: 



baryon acoustic peak measu rements (see, e.g 

2007 ; Parkinson et al]|2007 ). and galaxy cluster number counts (see, e.g.,|Voit 
20051 ) . may be used to constrain model parameters 



20071), including 



Doran et al. 



2005 



Younger et al. 



Since most observables depend on combinations of cosmological parameters rather then 
on just a single parameter, a single data set can not provide strong constraints. To get 
around this it is important to consider many different cosmological tests. This allows for 
consistency checks and might also allow for identification of systematic effects present in a 
particular data set. Combining data sets with constraints that are orthogonal to each other 
in parameter space results in significantly tighter constraints. 



In this p aper we use galaxy clu ster gas mass fraction versus redshift data (I Allen et al. 



2004 , also see Sasaki 



2005 



also see 



1996 



Penll 19971) and Hubble parameter versus redshift data (jSimon et al. 



Jimenez fc Loebll2002l ) to jointly constrain parameters of three different dark 



energy models. The first model we study is the cosmological constant dominated cold dark 
matter model (ACDM) with redshift-independent cosmological constant energy density pa- 
rameter f^A- We also consider the XCDM parametrization of dark energy, where dark energy 
is taken to be a fluid with an equation of state that relates pressure p x = ui x p x to the energy 
density p x , where uo x is a negative constant (this is only an approximate parametrization of 



1 Alternatively , it could be that general relativity needs to be modified on very large scales (see, e.g., 
Wang et al.ll2007t iMovahed et~al1l2007t lTsuiikawa|[iou7l ; lElizalde et aL1l2007h . 
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dark energy). Thirdly, we consider a slowly- rolling dark energy scalar field model (0CDM) 
in which the scal ar field <f> has potent i al energy density Vj^j) oc <fi~ a , where a is a nonneg- 



ative constant (IPeebles fc Ratralll988l : iRatra fc Peebles! Il988f ). For the </>CDM and XCDM 
cases we only consider spatially-flat spacetimes, while in the ACDM model spatial curvature 
is allowed to be nonzero. XCDM and 0CDM reduce to the time-independent dark energy 
ACDM model when u x = — 1 and a = 0, respectively. In this paper we jointly analyze both 
data sets and derive constraints on the nonrelativistic matter density parameter Q m and a 
parameter p that describes the DE. The parameter p is £7a for ACDM, u x for XCDM, and 
a for 0CDM. 

The galaxy cluster gas mass fraction versus reds hift data has been used to constrain pa - 
rameters of the ACDM, XCDM and 0CDM m odels (Ullen et al.ll2004 j: Ichen fc Ratral l2004f). 



These data provide tight constraints on fl m . iRapetti et al. 



used the galaxy cluster 



data in combination with CMB anisotropy and SNIa measurements to constrain dark energy 
evolution. For the XCDM model, assuming a time-independent equation of state, they set 
tight limits, u x = -1.05t° H, while more gen erally they fo u nd no significant evidence for 
evolution in the dark energy equation of state. I Wilson et al.l (120061 ) used these data in com- 
bination with SNIa data and found that the joint constraints were significantly tighter then 
those derived from either data set alone; the combined an alysis favored th e ACD M model 
but did not strongly rule out slowly-evolving dark energy. lAlcaniz &: Zhul (120051 ) used the 
galaxy cluster data and SNIa data (along with priors on the Hubble parameter and the bary- 
onic matter density) to jointly constrain brane world models. This data set has been used in 
conjunction with Fanaroff-Riley type lib radio galaxy angular size distance measurements to 
put an upper limit o n the amplitude of nq n-Riemannian terms during the late stages of the 
Universe's evolution (IPuetzfeld et al.ll2005l ). Galaxy cluster g as mass fraction data have also 



been used to constrain other dark energy models (see, e.g., IChang et al.l 120061 ; IZhao et al. 



200J). 



The H(z) data were used by lSamushia fc Ratral (120061 ) to constrain cosmological param- 
eters in the ACDM, XCDM and 0CDM models, but a comput ational error was made when 
cosmological parameter confidence contours were calculated. ISen &: Scherrerl (120071 ) used 
these data to constrain the evolution of an arbitrary dark energy component that satisfies 
the weak energy condition, in spatially-flat models. Th e H (z) data set has a ls o been used to 
const rain a number of interacting dark energy models ( Wei fc Zhang 2007al lbl; Izhang &: Zhu 
20071 ). In combination with CMB anisotr opy measuremen ts and SNIa data it has been 
used to constrain the Chaplygin gas mode l (IWu fc Yull2007l) as well a s cosmological models 
motivated by higher dimensional theories (ILazkoz fc Majerottol 120071 ). 



In this paper we present corrected cosmological parameter constraints for the H(z) 
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data. We also provide joint constraints on the ACDM, XCDM, and 0CDM models from the 
H(z) and galaxy cluster gas mass fraction versus redshift data. In Sec. 2 we outline our 
computational method. Results are presented and discussed in Sec. 3. 



2. Computation 



We use the I Allen et al.l (120041 ) measurements of gas mass fractions for 26 relaxed rich 
clusters in the redshift range 0.08 < z < 0.89. The cluster baryon mass is dominated by 
the gas. In relaxed rich clusters the baryon fraction should be independent of redshift. 
The cluster baryon fraction value depends on the angular diameter distance, so the correct 
cosmological parameter values place clusters at the right angular diam eter distance to ensure 
the redshift independence of the cluster baryon fraction. We follow IChen fc Ratral (120041 ) 
and compute the two dimensional likelihood function L G (Q m ,p) for each of the three DE 
models. When computing L (Q m ,p) we marginalize over the Gaussian uncertainties in the 
bias factor b, in the Hubble constant h ( in units of 100 km s~ 1 Mpc~ 1 ), and in the baryonic 
matter density parameter fib- Following lAllen et al.l (120041 ). we use b = 0.824 ± 0.089 (one 
standard deviation error) for the bias factor. To reflect the range of uncertainties, we use two 
sets of values for h and Q h h 2 . One set is Q h h 2 = 0.014 ±0.004 (one standard deviation error, 



Peebles fc Ratra 



Chen et al. 



2003] ) and h = 0.68 ± 0.04 (one standard deviation error, iGott et all 12001 
2003|). The other is from the WMAP t hree-year data, VL h h 2 = 0.0228 ± 0.0007 



and h = 0.73 ± 0.03 (one standard deviation errors, ISpergel et al.l 120071 ). 



The second data set we use are the nine lSimon et al.l (12005) measurements of the Hu bble 



parameter in the redshift range 0.09 < z < 1.75. Following iSamushia fc Ratral (120061 ) we 
compute a two dimensional likelihood function L H (Q m ,p) for each DE model. H(z) is not 
sensitive to the bias factor or baryonic matter density, but we still have to account for 
uncertainties in the Hubble constant. For the Hubble constant prior probability distribution 
function we use the same set of values as in the previous paragraph. 

To derive joint constraints, for each DE model we define the joint likelihood function 
L(Q m ,p) = L G (O m , p) L H (Q m , p) . From the joint likelihood function we compute 1, 2, and 3 
a confidence contours, as the contours that enclose 68, 95, and 99 % of the total probability. 



3. Discussion and Conclusion 



Figures [T] to [3] show cosmological parameter confidence contours for the ACDM, XCDM 
and 0CDM models for the two sets of Vt^h 2 and h priors. 
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Figure [T] shows constraints on the ACDM model. The galaxy cluster gas mass fraction 
data place a good constraint on Q m (< 0.35 at 3 a), while the H(z) data constrain a linear 
combination of fi m and The joint likelihood functions peak near spatially-flat models. 

Figure [2] shows the constraints for the XCDM parametrization. The joint constraints 
favor the region of parameter space near the u x = — 1 line which corresponds to spatially-flat 
ACDM models. 

Figure [3] is for the 0CDM model. The joint likelihoods peak on the a = line which 
corresponds to the spatially-flat ACDM model. However, values of a as high as 4 or 5 are 
allowed at 3 a. 

The galaxy cluster gas mass fraction data is more restrictive than the H(z) data. When 
they are combined the H(z) data shifts the constraints to slightly higher values of Q m than 
for the galaxy cluster gas mass fraction data set alone. A spatially-flat cosmological model 
with a cosmological constant term with Q\ ~ 0.7 is a good fit to the joint data in all six 
cases consider e d her e. This is consistent w ith re sults based on othe r measurements, see, e.g., 
Rapetti et all j2005[ ). Iwilson et all J2006h . and bavis etaJ J2007h . 



Hubble parameter versus redshift data is expected to increase by an order of magnitude 
in the next few years. In combination with new galaxy cluster gas mass fraction, SNIa, and 
CMB measurements, this will significantly better constrain dark energy models. 

We thank R.Lazkoz for helpful discussions. We acknowledge support from DOE grant 
DE-FG03-99EP41093, INTAS grant 061000017-9258 and NASA ATP grant NAG5-12101. 
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Fig. 1. — 1, 2, and 3 a confidence level contours for the ACDM model. Dashed lines denote 
constraints from Hubble parameter versus redshift data, while solid lines show the joint 
constraints (the crosses indicate the maximum likelihood points). The diagonal dotted line 
corresponds to spatially-flat ACDM models. Thick lines correspond to the h = 0.73 ± 0.03 
and Vt h h 2 = 0.022 ± 0.0007 priors (maximum likelihood is at fl m = 0.26 and fi A = 0.85), 
while thin lines are for h = 0.68 ± 0.04 and Vt^h 2 = 0.014 ± 0.04 (maximum likelihood is at 
Q m = 0.18 and tt A = 0.70). 
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Fig. 2. — 1, 2, and 3 a confidence level contours for the XCDM parametrization. Dashed 
lines denote constraints from Hubble parameter versus redshift data, while solid lines show 
the joint constraints (the crosses indicate the maximum likelihood points). The dotted 
horizontal line corresponds to spatially-flat ACDM models. Thick lines correspond to the 
h = 0.73 ± 0.03 and Vt h h 2 = 0.022 ± 0.0007 priors (maximum likelihood is at fi m = 0.26 and 
cj x = -1.2), while thin lines are for h = 0.68 ± 0.04 and Vt h h 2 = 0.014 ± 0.04 (maximum 
likelihood is at f2 m = 0.20 and uj x = —0.98). 
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Fig. 3. — 1, 2, and 3 a confidence level contours for the 0CDM model. Dashed lines denote 
constraints from Hubble parameter versus redshift data, while solid lines show the joint 
constraints (the crosses on the horizontal axis indicate the maximum likelihood points). The 
horizontal a = axis corresponds to spatially-flat ACDM models. Thick lines correspond to 
the h = 0.73 ± 0.03 and Vt h h 2 = 0.022 ± 0.0007 priors (maximum likelihood is at fi m = 0.26 
and a = 0), while thin lines are for h = 0.68 ± 0.04 and Vt^h 2 = 0.014 ± 0.04 (maximum 
likelihood is at f2 m = 0.20 and a = 0). 



